Universality in p-spin glasses with correlated disorder 
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We introduce a new method, based on the recently developed random tensor theory, to study 
the p-spin glass model with non-Gaussian, correlated disorder. Using a suitable generalization of 
Gurau's theorem on the universality of the large N limit of the p-unitary ensemble of random 
tensors, we exhibit an infinite family of such non-Gaussian distributions which leads to same low 
temperature phase as the Gaussian distribution. While this result is easy to show (and well known) 
for uncorrelated disorder, its robustness with respect to strong quenched correlations is surprising. 
We show in detail how the critical temperature is renormalized by these correlations. We close with 
a speculation on possible applications of random tensor theory to finite-range spin glass models. 



Introduction. It is well known that the phase diagram 
of the p-spin glass model [H, 0] does not depend on the de- 
tails of the disorder distributions, in the following sense: 
CZ2 ■ if Ji t ...i denotes a set of independent and identically dis- 



tributed set of p-valent coupling between sites i\ 



non-quadratic potential V(J» I ...i 
bution 



in the coupling distri- 
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is irrelevant in the thermodynamic limit. That a similar 
result would hold for correlated disorder distributions, 
with terms such as 
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in the potential, is much less obvious. In fact, to our 
knowledge, no analytic framework to deal with such cor- 
related, non-Gaussian disorder has been reported so far. 
Since disorder correlations are to be expected in actual 
physical systems, understanding their effect is an impor- 
tant problem. 

In this Letter we exhibit an infinite class of non- 
Gaussian terms of the kind (J2]) such that (i) the ther- 
modynamic limit N — > oo is exactly soluble, and (ii) the 
spin glass phase has the same structure as with uncorre- 
lated disorder, except for a renormalization of the critical 
temperature. This provides the first general result on spin 
glasses with strongly correlated disorder. 

Our approach is based on new results in random tensor 
theory. As natural generalizations of random matrices, 
random tensors have recently been showed to possess a 
large iV limit Q dominated by few, well-identified "mel- 
onic" graphs (the tensor equivalent of 't Hooft's planar 
graphs in matrix theory Q). Furthermore, the melonic 



family can actually be resummed exactly, and turns out 
to exhibit interesting critical and multicritical behavior 
These results have not been applied to spin-glass 
problems previously, and our hope is to convey that ran- 
dom tensors are potentially as powerful tools for spin 
glass theory as random matrices [6[. 

From the perspective of random tensor theory, the 
quenched couplings of spin glasses with p-spin interac- 
tions are non-Gaussian rank-p random tensors. The be- 
havior of such tensors in the large N limit has been inves- 
tigated in 0, [H , with a striking conclusion: in a suitable 
ensemble with p-unitary symmetry (more details in the 
text), this limit is universally Gaussian. This means that, 
in this ensemble, in the large TV limit the sole effect of the 
self-interactions of large tensors is to dress the propaga- 
tor. Here, we show how this result can be generalized to 
include interactions between tensors and spin variables, 
and thus obtain the aforementioned universality result. 

This Letter is organized as follows. We first recall the 
Hamiltonian for p-spin models and insist on the need for 
a correlated disorder. Then, we recall the relevant prop- 
erties of large random tensors in the p-unitary ensemble. 
This enables us to show how non-Gaussian, correlated, 
quenched variables can be integrated exactly in the large 
N (thermodynamic) limit, yielding our universality theo- 
rem. We conclude with a few words on the possible rele- 
vance of tensor techniques for short-range p-spin glasses. 
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We consider a p-spin Hamilto- 



Hj(S) 
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where Ji 1 ...i is a com ple^fl tensor describing the cou- 
plings and S = (Si)i<i<N is a set of real spins@ with 
lattice index i, weighted by a (normalized) probability 
measure dQ(S) such that 



/ 
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This includes in particular Ising jl| and spherical Q spins. 

When the couplings are Gaussianly distributed, such p- 
spin glass models are well-known to exhibit replica sym- 
metry breaking in the low temperature phase 0, [H , 12 1 
and have a dynamical transition at a higher temperature 
where a large number of mctastable states (growing ex- 
ponentially with N) dominates the free energy landscape 
13l . Il4| ; their relevance is conjectured to extend to struc- 
tural glasses [lH . These results extend easily to the case 
of independent and identically distributed (i.i.d.) cou- 
plings: all terms of higher order than 2 in Ji x ...i and 
Jii...i p in the measure on (J, J) are irrelevant in the ther- 
modynamic limit ./V — > oo. 

In this Letter we aim to study a family of conflated 
non-Gaussian measures on the disorder. Physically, ran- 
domness of the couplings comes from randomness of the 
positions of the spins, and in general we should not ex- 
pect the couplings between different sets of p spins to be 
independent (for instance due to the geometric relations 
between the positions of the spins). One should therefore 
perturb the Gaussian distribution (with covariance a 2 ) on 
the couplings with a polynomial V{J, J). The quenched 
free energy is given by 



[F{J,J)\ 
with 
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and J • J is shorthand for Y] 

The evaluation of [F] for a generic potential V is of 
course a completely open problem. However, the theory 
of large random tensors provides an exact calculation for 
an infinite family of potentials satisfying a particular kind 
of invariancc. 



1 The use of complex rather than real tensors is motivated by 
purely technical convenience and does not change the physics 
in any way. 

2 It is also possible to include Potts or vector spins coupled accord- 
ing to some fixed multi-linear map. 
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FIG. 1. Some p-bubbles at p = 3. Up to color re-labeling 
there is a single bubble with 4 vertices (on the left), whose in- 
variant is ^{ij.j,} Jiii2*3 J'ij233 J313233 J iii^z • But there exist 
different bubbles with 6 vertices (the three on the right). 



Large random tensors. We now review the relevant 
properties of large random tensors discovered in 0, 0, 0] ■ 
The first obvious observation is that, unlike symmet- 
ric/hermitian matrices, tensors cannot be diagonalizcd. 
Hence, a key concept in random matrix theory, the eigen- 
value distribution, does not carry over to the higher-rank 
case. It turns out however that this fact does not preclude 
the development of random tensor theory, which in fact 
relies on the identification of an ensemble with suitable 
symmetry properties. 

One such ensemble of tensors — indeed the only one 
identified so far — is the p-unitary ensemble, defined as 
follows. Consider a rank-p tensor in N complex dimen- 
sions J, with components Ji t ...i in a fixed basis, and for 
each set of p unitary matrices to U^ p \ define 



J' 
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Then we say that a function V(J, J) of J and its complex 
conjugate J is a p-unitary invariant i^| 



V(J',J') = V(J,J). 



(8) 



The set of p-unitary invariants is conveniently 
parametrized by p-bubbles B, that is p-valcnt bi- 
partite connected graphs with edges colored by numbers 
between 1 and p, such that each "color" is incident 
exactly once to each vertex, sec Fig. [TJ A bubble 
represents an invariant denoted tr^(J, J), by associating 
a tensor J to each "white" vertex of B and a conjugate 
J to each "black" vertex, and contracting their fc-th 
indices along the edges colored by k. By the fundamental 
theorem of classical invariants of U(N) (see for instance 
0]), a general p-unitary invariant can be expanded as 



V(J, J) = J2t B tT B (J, J), 



(9) 



3 The corresponding symmetry group is known as the external ten- 
sor product of p copies of U(N). 



where ts are coupling constants. 

For a given invariant potential V(J,J), we define the 
average of f(J, J) over J by 
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The Fcynman diagrammatic expansion of these quanti- 
ties involves (p + l)-colored bipartite graphs, made of p- 
bubblcs connected together via extra lines with color "0" 
incident on each vertex and corresponding to the propa- 
gator a 2 in ([TO]) . 

The following results concerning the large N limit of 
(fT0|) have been proved: 

• The Feynman expansion is dominated in the large 
N limit by a simple class of graphs, called melonic 
graphs, which generalize 't Hooft's planar graphs 
[3j. Intuitively, a (p + l)-colored graph is melonic 
if it can be built by recursive insertions on any line 
of two vertices connected together by p lines, as in 
Fig.H 

• The large N limit is Gaussian, in the sense that up 
to subleading corrections in 1/N, 



tr B (J, J) 



NG 



\B\/2 



(11) 



where \B\ is the number of vertices of the bubble 
B and G2 = [J ■ J)/N is the dressed propagator 
depending on the potential V 

• The following Schwinger-Dyson equation holds in 
the N -> 00 limit [§] 



[J- J] , ^ \B\ [tT B (J,J)] 



a 2 N 



N 



(12) 



The first result implies that all non-melonic bubbles B in 
the potential drop out in the large N limit, and therefore 
we can restrict the sum in ((9|) to melonic bubbles (hence 
hereafter B will always denote a melonic bubble). In Fig. 
[TJ all bubbles are melonic except the non-planar one on 
the right. 

The second result has been coined the universality 
property of the p-unitary ensemble of random tensors, 
and can be seen as a non-trivial generalization of the 
central limit theorem. Its origin is that there is only one 
way to dress a melonic bubble B with propagators in a 
melonic way, which happens to correspond to Gaussian 
contractions. This feature is specific to tensors and does 




FIG. 2. A patch of a melonic graph, recursively built by in- 
serting on any line a pair of vertices, a black and a white, 
connected together via p lines (here p — 3) . 



not hold for random matrices. In a way, this Letter can be 
read as the physics counterpart of this surprising mathe- 
matical result. We refer the reader to the review and 
to the original papers for more details on random tensor 
theory. 

Universality in the couplings. Let us now come back 
to spin glasses. Followin g th e standard recipe to com- 
pute quenched quantities [17| . we consider the averaged 
replicated partition function 



/n 
~[dn{s a ) e - 0H ^ 



as a }) 



(13) 



where a is the replica index and the effective Hamiltonian 
is defined by 



-/3H eB ({S°-}) 



(14) 



In diagrammatic language, H c g is given by the sum over 
all connected (p+ l)-colored bipartite graphs (henceforth 
"graph") with spins Sf on the external legs. Denoting 
k the order of the effective coupling between replicas 
a\ ■ ■ ■ a/., this can be pictured as 




, a 2 



Gl 



ai...o fc <2{: 



(15) 

Here the solid line is the J-propagator (tensor lines with 
color 0), and the p dashed line emerging from each ex- 
ternal leg represents the external spin variables S'° fc . The 
blob Gk is the large- N tensor connected fc-point function, 
i.e. the sum over all connected melonic graphs with k ex- 
ternal (solid) legs. For each graph contributing to the 
blob amplitude, the site indices ii of the spins are con- 
tracted along "broken faces", i.e. connected paths with 
alternating color 1 < c < p and from one external 
dashed leg to another through the graph. 
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Let us now show that k = 2 terms dominates in the 
large TV limit. Observe that powers of TV in H c s({S a }) 
have three sources: the tensor propagators, the bubble 
interactions tr#(,7, J), and the sums over site indices i. 
The first two contributions are those of a melonic graph 
with k cut lines of color whose scaling has been was 
found in the appendix of to be p — (p — l)k — p, where 
p is a positive number independent of A:. As for the spin 
contribution, from (j4]) we see that it gives at most a factor 
of TV per broken face, and there are at most pk/2 of them. 
This gives for the scaling degree cj(k) in TV of the order-A: 
term of (fT5]l 



w(fc)<p-(|-l)fc. 



(16) 



We conclude that, indeed, only k = 2 terms are relevant 
in the large TV limit. Thus, at leading order (fTijj) reduces 
to 



/3iJ off ({5 a }) = 



a.b 




(17) 



To complete our evaluation of the effective Hamil- 
tonian, we must compute the 2-point function 
(G2)ii-i p ;h-j p = [Jii-i p Jji-j v ] 01 tnc tensor. Its 
scaling with TV is TV~( p_1 ). Its tensorial structure is 

n 



1=1 u U-]l 



which identifies by pairs the lattice sites 
between the replicas a and b. Finally its amplitude, 
simply denoted G 2 , is found by inserting the universality 
property (fTTj) into the Schwinger-Dyson equation (fl2"j) . 
yielding 



rn>2 V SSB™ 



t B 



mG? = 1, 



(18) 



in which B m denotes the set of melonic bubbles with 2m 
vertices. The leading-order connected 2-point function is 
the solution of this polynomial equation, and depends 
on the whole set of coupling constants ts- For example, 
for a potential with a single 4-vertex bubble (see Fig. [T]) 
with coupling constant t, equation (|18[) becomes 

2a 2 tG 2 (t) 2 + G 2 (t) =ct 2 , (19) 

hence, picking the solution with G2(0) = a 2 , 



G 2 (t) = 



VI + 8<jH - 1 



(20) 



This is a smoothly decreasing function of t > — l/8tr 4 . 
Summarizing, we have proved that 



fi 2 G 2 \ -> \ - 

a,b i\---ir> /— 1 



b 



(21) 



which is the usual p-spin replica Hamiltonian (ill . [lij , ex- 
cept for the variance a 2 which is replaced by G2 (which 
as we saw can be computed exactly for a given tensor 
quenched potential V). This is the content of our univer- 
sality theorem, the main result of this Letter. It shows 
that the higher order terms in the quenched distribution 
change the critical temperature, but not the structure of 
the low temperature phase. 

Conclusion and outlook. We have introduced large 
random tensors as a new tool for spin glass theory. Using 
the peculiar scaling behavior of tensors in the p-unitary 
ensemble, we have identified an infinite universality class 
of infinite-range p-spin glasses with non-Gaussian corre- 
lated quenched distributions. To our knowledge, this is 
the first universality theorem of spin glass theory with 
this level of generality. 

We close with a more prospective remark. Just like 
their Shcrrington-Kirkpatrick relatives, the p-spin inter- 
actions in ([3]) have infinite range, and for this reason p- 
spin glass models are judged (at least partially) unphys- 
ical. We expect however that random tensor techniques 
should be applicable to finite-range models too. Indeed, 
from the random tensor perspective, a finite-range spin 
glass model is one for which the J-propagator is non- 
trivial, and in particular depends on the tensor indices of 
J. A typical example of interest here would be 



Ex 



<l<k<p 



(ll-Jk) 2 + 1' 



(22) 



Or analytic, if V has infinitely many bubble terms. 



which goes to zero when the lattice sites are far away. 
Such tensor models have already been considered in the 
context of quantum gravity [18(, where they have been 
called tensor field theories (TFT). The key difference be- 
tween TFT and the simple tensor models considered in 
this Letter is the appearance of a renormalization flow. 
The first renormalizable TFT has been identified in fl9j ] . 
and developments are fast in this area. We expect that 
these new techniques will prove useful in the difficult field 
of finite-range spin glass theory. 
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